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Scope and Purpose—A central planning problem in both the public and private sectors is determining the 
best locations for a set of facilities. Each facility is assumed to house a number of servers who provide 
some service to a spatially dispersed set of customers. In many systems, customers are not scheduled— 
they arrive randomly in time—and the amount of service they require is also an uncertain quantity. In such 
cases, facility congestion can occur in the sense that all servers associated with one or more facilities become 
busy. While a considerable amount of effort has been devoted to optimal facility location, most of it has ignored 
this type of congestion, assuming essentially unlimited service capacity at each facility. In this paper we 
explicitly take facility congestion into account by allowing a customer to be serviced by servers from other than 
the closest facility—should the closest be saturated. In fact, depending on the pattern of facility congestion. 
any given customer may be serviced by a server from any one of the facilities. Given the model assumptions, 
the objective is to locate the facilities on the associated transportation network so as to minimize mean travel 
time to a random service request. This work begins to tie together previously disparate work in (deterministic) 
location theory and spatially distributed queues, In particular, the results are applicable to the hypercube 
queueing model, first reported seven years ago in this journal. 


Abstract—The median problem has been generalized to include queueing-like congestion of facilities (which 

are assumed to have finite numbers of servers). In one statement of the generalizations, a closest available 
server is assumed to handle each service request. More general server assignment policies are admissable. 

The objective is to minimize the steady state expected travel time associated with a random service request. 

It is shown that, under suitable conditions, at least one set of optimal locations exists solely on the nodes of 

the network. It is also shown that this result has a direct relationship to the hypercube queueing model. 


INTRODUCTION 

The problem of determining the locations of a set of facilities on a network that minimize the 
expected travel time to or from the facilities, for the population of their users, is one of the 
classic problems in location theory. This problem, known in the literature as the median 
problem, has been studied very thoroughly in the last two decades. The basic theoretical results 
in this area are due to Hakimi[5, 6]. Subsequently, Goldman[4], Hakimi and Naheshwari[7], 
Levy[11] and Wendell and Hurter[13] have extended and generalized Hakimi’s results. 

When there are Q facilities to be located on an undirected network G, the median problem is 
to find a set Z* =(Z¥, Z3,..., Z§) of Q points on G such that 


S hd(Z*, = hdlZ, (1) 
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for all Q-tuples Z of points on G, where fh; is the fraction of demand that is generated at node 
j 1 4j=1), n is the number of nodes and d(Z, j) is the shortest distance between node j and 
the closest point in Z. In [5] Hakimi proved that at least one such set Z* exists solely on the 
nodes of the network, so that the search for an optimal solution can be confined to the finite 
node-set of G. 

For this standard median problem to be applicable, four main assumptions should be 
satisfied: 

(1) Travel in the given area is restricted to take place solely along the links of the network. 

(2) Requests for service can occur only at the nodes of the network. 

(3) When the number of facilities is greater than one, a service request from a particular 
location is always handled by a server at a closest facility. 

(4) There is always an available (free) server at the selected (closest) facility. 

The traveling associated with a service request could require the “customer” (requester of 
service) to travel to a nearest facility or a server at a nearest facility to travel to the customer. 
The former, ‘“‘customer-to-server” type system, includes outpatient clinics, “little city halls”, 
libraries and even hamburger havens. The latter, “‘server-to-customer” type system, includes 
emergency services (e.g. police, fire, ambulance, emergency repair), special-order delivery 
services, and certain home visitation medical services. In our work, we use the term “travel 
time associated with a service request” to mean either the customer-to-server or server-to- 
customer travel time. 

The type of systems we consider are characterized by stochastically generated requests for 
service (in time and space) and by nondeterministic service times for meeting these requests. [A 
service time is comprised of travel time plus on-scene time.] In such an environment, it is quite 
possible that all servers at a nearest facility will be busy, thereby yielding a congested network 
in which queues could form. Thus, assumption (4) above often does not hold. For these 
systems, equation (1) expresses only an unnatural problem: finding a set of points so as to 
minimize the expected travel time for a random service request arising at one of those special 
times when servers are available at all the facilities. In contrast, it is the purpose of this work to 
incorporate in the median problem the possibility that all servers at any subset of the Q 
facilities can be busy. 

The objective function in this median problem with congestion is to minimize the expected 
response time associated with a random service request, where response time is the sum of 
travel time and queueing delay. The averaging must be carried out using the equilibrium state 
probabilities of the system, defining “‘states” according to the status of each of the facilities—at 
least one server available at the facility or all servers busy. To avoid queue formation whenever 
possible, we assume that a most preferred available server is immediately dispatched to any 
service request. A request is placed in queue only if all servers are simultaneously busy; such a 
queue is assumed to be depleted in a first-in, first-out (FIFO) manner. Usually, server 
preferences are dependent solely on geographical proximity, but more general server pref- 
erence policies are allowed. 

In a recent paper[1] we solved this problem in its full generality for the case of one facility 
staffed with one server. Using results from the theory of M/G/1 queues, we developed an 
exact procedure for finding the optimal facility location, which was found to be sometimes ona 
node and oftentimes on a link. This result, which contrasts sharply with the nodal solution 
result of Hakimi and others, was due to nonlinearities in the mean queueing delay expression. 

Our purpose in this paper is to extend the 1-server results to the case of multiple servers 
(and multiple facilities). However, due to the analytical intractability of M/G/n queues 
(m =number of servers = Q), we have only been able to develop results for a special case, 
namely the case in which on-scene service time is much greater than travel time—so that state 
probabilities depend only on three sets of input parameters: nodal service request rates, 
on-scene service times, and server preferences. Changes of facility location occurring with all 
three parameter sets fixed do not change the state probabilities. However, this does not 
preclude dependence of system state probabilities on facility location: any change of facility 
location that creates a change in server preferences will, in general, change the state prob- 
abilities. With these restrictions, we prove a nodal result analogous to Hakimi’s. However, it 
should be clear from the 1-server case that if the service time assumption is not valid in a 
particular setting, then optimal server locations are not in general at nodes. 
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The analysis, given the aforementioned restrictions, ties together previously disparate 
research efforts on network analysis and on spatial queueing analysis. In particular, we show 
that the hypercube model[9, 10] and the algorithm of Jarvis[8] for determining optimum 
locations can be useful to solve the median problem with congestion for specific situations. In 
addition this work indicates that the basic hypercube model does not suffer a loss of generality 
by considering only nodes (or atoms in the context of the hypercube model) for the locations of 
the service units. 


NOTATIONS AND ASSUMPTIONS 


Let G(N, L) be an undirected network where N is the set of nodes with |N|=n, and L is 
the set of the links. Let X be the set of all possible locations of Q facilities (Q > 1), on the 
network G, i.e. 


Kee oH Giessole REE “hehe Qk: 


Given any Q locations Xg = (i,,.. sig) € Xe let i denote the event that the facility at i, is 
not staffed with an available server (the facility is busy) and j, that the facility at i, does 
have an available server. Therefore, for any Xp € Xq there are 2° combinations (states) of 
finding the network at any time, according to the status of the Q facilities. Let Yx;q) be the set 
of all states for X9EXo and let yx(q (or for convenience yg) be a generic element of ¥x(Q) 

We assume that server assignment occurs according to a fixed preference procedure. That is, 
for each demand point in the network there is a list of facilities that specifies the ordering of 
preferences for the assignment of servers (i.e. first preference for servers from facility k,, 
second preference for servers from facility k,, etc.). A most preferred available server is 
always assigned to a customer.t The goal of the optimization to be stated below is to minimize 
expected system travel time under a given fixed preference procedure. The fixed preference 
itself need not be determined solely by relative travel times, but can include individual 
characteristics of servers (e.g. bilingualness) and special needs of customers at the nodes. 

Let ¢(i, j) be the travel time on link (i, j), (i, J) € L, and let d(yg, j) be the (minimum distance 
path) travel time associated with a most preferred available server to node j, when the system is 
in state yo. 

As in the standard median problem we assume that service requests are generated on the 
nodes of the network. However, in addition, we assume that service requests occur according 
to a homogeneous Poisson process, with each request requiring a service time whose dis- 
tribution is general but not dependent on its location or the location and identity of the server or 
the history of the system. This implies that on-scene service times are assumed to be much 
greater than travel times so that variations in total service times which are due solely to 
variations in travel times among potential servers are ignored. (This assumption is reasonable 
for systems whose on-scene service times are at least an order of magnitude greater than travel 
times.) 

Finally, we require that travel time is uniform over a link, i.e. the travel time over a fraction 
6 of some link (p, q) is 6t(p, q). This assumption is not restrictive since the links and nodes can 
be defined in such a way that this assumption holds to a specified degree of accuracy. 


MODEL FORMULATION AND ANALYSIS 


We will consider the steady state behavior of the system. For any possible set of locations 
Xq € XQ, let P(yg) be the steady state probability that the network is in state yg © Yucq). (We 
assume that the appropriate ergodicity conditions apply so that a unique steady state dis- 
tribution exists.). Let yg’ be the state in which all the Q facilities are busy (ie. yg’ = 
(i, b,..., @) in our notation). 


tWhen preferetices depend directly on travel times, such a myopic strategy (through very reasonable) is not always 
optimal in the sense of minimizing time-average mean travel time. An optimal policy occasionally requires assignment of other 
than the most preferred available server [2], in order to leave the system ina state which best anticipates future service requests. 
We do not consider such strategies in our formulation of the median problem with congestion. 
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Conditioned on any state yg © Yx<q)— {Yq’}, the expression 
> hjd(yo, ) 
t 


is the expected travel time associated with a random service request. 
Suppose now that the network is in state yo’. We will consider three alternative policies 


regarding this state: 


(a) Service requests that occur while all the service units are busy, are handled by a back-up 
service system (zero line-capacity case). Let R be the travel time cost of utilizing this special 
reserve service system. 

(b) Service requests that arrive while all the facilities are busy enter an infinite capacity 
queue that is depleted in a first-in, first-out (FIFO) manner; upon completion of service, the 
server is either assigned to the next request waiting in queue, or returns immediately home if 
none is waiting. Therefore, if d(k, j)= travel time between nodes k and j, 


Dy > hyhyd(k, j) 


is the expected travel time of a random service request given that the network is in state yo’. 

(c) Requires negative exponential service times. As in the previous case, service requests 
that arrive while all the facilities are busy enter a FIFO queue with infinite capacity, but now 
upon completion of service, the server always first returns to his/her home location. In this 
case, assuming negative exponential service times, the conditional expected travel time of a 
random request is 


SSO nai, 


given that the network is in state yg’, where n(k) = number of servers stationed at facility k, 
located at i,, and 7=Z2, 7(k) is the total number of servers. 

The appropriateness of any particular assumption depends, of course, on the system being 
modeled. Assumption (a) often applies to ambulance systems, in which emergency requests cannot 
be queued. Assumption (b) applies frequently to police vehicles that may be dispatched 
back-to-back to successive service requests. Assumption (c) applies to some ambulance and fire 
services, but suffers from the assumption of negative exponential service times. The median 
problem with congestion is now stated: 


min F(Xg) (2) 
XQEXQ 
with 
F(X)= SY , Pla) D hiya i) + Pye") & hiCU) 
YQEYx(Q)-{¥Q"} j=l a 
where 


CU) is R or > h,d(k, j) or > -_ (ix, j) 


according respectively to (a), (b) or (c) above. 

Obviously, the standard median problem is a special case of (2) arising when P(yg)=0, 
V yo#(i,...,/g)—the state where all the facilities have one or more units available and when 
d(yg, j) is determined solely by geographic proximity (i.e. minimizing travel time). In contrast, 
the weights P(yg) in (2) represent the fractions of time that the network is in each of the 22 
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possible states. Therefore, as noted before, we take into account that any subset of facilities 
can become depleted of servers. 


Now, given the assumptions stated before, the following theorem can be proved. 


Theorem 

For a given fixed preference server assignment procedure, at least one set of optimal 
solutions to (2) exists on the nodes of the network. 

Proof. Let Xg=(i,, iz,..., ig) be any feasible solution to (2), and let P(yg), Vyg€ 
Y xq) be the corresponding steady state probabilities. Suppose that is is an interior point on the 
link (p, q). Then by the uniform speed assumption 


t(p, ig) = r) 


t(D, 4) 0<6<1. (3) 


The following proof is for the case C(j) = 22, (n(k)/) d(i,, j) in (2). The proofs for the other 
two cases are very similar and even slightly easier. Let Y;,C Yx:q)—{Yg’} be the set of all 
states in which the facility located at i, is available (i.e. has at least one available server). Then 
we can write F(X) as: 


FIXe)= Pla) hdiveri+ Pose’) { SMO hacia] +4 (4) 


YQE Vig 


where the term A includes all server assignments that exclude the facility located at ig, i.e. 


yaE¥xQ)- Yig-lyQ?} 


A= SPS hdvod+POe[S & Madina). 
KES 


Let N,,(is) be the set of all nodes that would be eile to a server from the facility at ig, 


when the network is in state yo € Y;,, and let N, glis) = N — Nyo(is)- Then we can rewrite (4) as: 


FUXa)= 5 POo)|,_F, hain ]+% POS hdiod] rar B — 


where the term B corresponds to non-queued assignment of servers not located at is, even 
when the facility located at i, is available, i.e. 


B= P(yg) SY  hjd(yvq/). (7) 


49E Nig JEN, lis) 


Recalling that is is assumed to be an interior point on the link (p, q), let N,o(is, p) C Nyg(is) 
be the set of all nodes that belong to the set N,,(is) and which communicate most efficiently 
with the facility at is via p, and let Nyo(is, g)=Nyo(is)—Nyo(is, p). (The term “communicate” 
implies minimal travel time.) If a node communicates equally ara with is via nodes p or q 
for some yo, we can include that node in either N,_(is, p) or N, (is, q), but not in both. 

Let N(is, p) be the set of all nodes which communicate mat “efficiently with the facility at is 
via node p and let N(is, q) = N — N(is, p). 

Therefore, we can write (6) as 


FIX)= B Povod{ Dill p)+ is) 


YQE Vis vals. p) 


+S, ,,hlddi a+ Hq isn} + 22 Pog Shlain) +e. i 


JEN, yolis. q) 


+ 3 Ahldi,a)+0g, ih} + A+B, 8) 


ieNUs, @) 
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Using (3) and rearranging terms we get 


FiXo)= a 10.0 B Poo) 3 htt Poe), By] 


IENyolis. P) 


+(1-8) [t. 9) (5 Ply) XS +12) py 0) A) | 


JENyolis. q) jEN(is, 
+A+B+C 9) 


where the term C corresponds to “fixed components” of travel time relating to the link (p, q), 
1.€. 


C= ¥ Pood{ 3D bdlir+ > i) 


YQE Vig FENyolis, P) JENyolis. 4 
Dry) FY ndir+ DY hada], a 
n jéEn(is, p) jENs, a) 


Since the steady state probabilities are independent of 6 (refer to the assumption from 
previous section), F(X) is a concave function of @. This implies that the minimum of F(XQ) as 
a function of @ occurs at an extreme point, either @=0 or 1, corresponding to locations at node 
p or q, respectively. Clearly, the node p is optimal if the coefficient of @ in (9) is larger than the 
coefficient of (1 — 6); otherwise g is optimal or a tie exists, in which case either is optimal. We 
showed so far that by starting with any Xg which is not at all nodes, there is an XQ! with 
F(Xq')<F(Xg) which agrees with Xg except that one non-node point is of Xg has been 
shifted to one of the two nodes adjacent to it. This argument can be repeated to show that there 
is an all-node XQ with F(X’) < F(X). Thus, the search for an optimum can be restricted to 
nodes; since there are only finitely many sets of Q nodes, the minimum must be attained 
(minimum instead of merely infimum). 


THE MEDIAN PROBLEM WITH CONGESTION AND THE HYPERCUBE MODEL 


“The hypercube model” is a spatially distributed queueing model developed by Larson[9] to 
investigate analytically the performance of urban emergency services. The model assumes a 
geographical region R that is divided into n geographic areas or atoms. The fraction of demand 
associated with each atom j is A,(27_, h; = 1) and the travel time from atom i to atom j is d(i, j). 
Service requests over the entire region are generated in a Poisson manner at a rate A and at 
each atom j independently in a Poisson manner with rate A,(2; A; = A). 

There are a total of Q service units (i.e. servers) to respond to the requests for service, 
located at atoms i,, i,,..., ig. For Markov analysis, the service time for each unit n is assumed 
to be exponential with mean y,,'. Recent research has shown that the assumption of exponen- 
tiality of the service time does not markedly affect the predictive accuracy of the model [8]. 

States of the system are defined by the status of each service unit: busy or available. The 
model allows a zero line-capacity (implying the existence of special reserve units), or an infinite : 
capacity queue. Given a dispatching policy, all the 22 steady state probabilities of the system 
can be obtained by solving 22 detailed balance equations [9]. In [10], Larson used a server 
sampling scheme adapted from the M/M/Q model to obtain fast approximate solutions for the 
required probabilities. 

For a given set of single server locations at atoms i,,..., ig the hypercube model computes 
several performances measures. Among them, the most imporant one is the region-wide mean 
travel time, defined as 


> > iyi Aix i) + P (all units are busy) 3 A, Cli) (11) 
7=1 k=1 id 


where p,.;=fraction of all nonqueued dispatches that send the unit from atom i, to atom j; 
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k=1,...,Q;j=1,...,2; C(=average travel time to node j arising from dispatches from queued 
service requests (infinite capacity case) or from service requests handled by a back-up service 
system. The p,j’s represent the response patterns of units. 

In [8] Jarvis developed an algorithim to find a set of “optimum” locations in the framework 
of the hypercube model where locations are constrained to atoms. The key idea behind the 
Jarvis algorithm is to optimally locate the servers (facilities) for a given response pattern and 
then, given a new set of locations, to reassess the response patterns to determine if a new set 
of dispatch preferences (and thus response patterns) could improve system performance 
further. This alternating iterative procedure is analogous to the “locate~allocate” scheme often 
used in deterministic location theory [12]. 

Jarvis’ algorithm for the zero line-capacity case works as follows: 


(1) Initialize: Specify initial unit locations for units 1,2,...,Q, corresponding to atoms 
i, b,..-, ig. 

(2) Allocate: Solve the hypercube model to obtain p,;,k=1,...,Q;f=1,...,m. 

(3) Locate: Solve the following L.P problem: 


min > SD Ptr, b) Cv, b) 
=loe=1 
st. ¥ P(v,k)=1 k=1,...,0 
v= 
P(v, ky) 20 v=1,...n;k=1,...,Q 


where the decision variable P(v, k) is the probability that server k is at node » when available, 
v=l,...,n; k=1,...,Q; and Clo, k)=22 p,;d(v,)), v=1,...,0; k=1,...,Q The LP 
has an optimal solution with P(», k)= 1 for some v and every k. 

(4) Test for convergence: If the new Q locations are identical to the old set of Q locations, 
stop. Otherwise, go to Step 2 with new unit locations i,,..., ig, and reallocate. 

Whenever the algorithm terminates, at least a local optimal solution is ensured. By taking 
several different initial sets of locations, the chances of getting closer to the optimal global 
solution are improved. 

The hypercube model can be applied in our congested median network context. The network 
G can represent the geographical region R, the nodes of the network being the atoms, and the 
links being the major streets connecting the atoms. It requires only some algebraic manipula- 
tions to demonstrate that if we take any Q points in the network to be the set of server 
locations, then F(Xg)—the cost function of our problem (2}—turns out to be identical to the 
mean region-wide travel time of the hypercube model [9]. 

The conclusion of this discussion is that since the assumptions of our Theorem hold for the 
hypercube model (subject to our discussion of service times), neither the hypercube model nor 
Jarvis’ algorithm suffer a loss of generality by considering only locations on the atoms. In 
addition Jarvis’ algorithm can be applied to the median problem with congestion whenever the 
hypercube model’s assumptions are accepted. This result ties together two very diwerent 
approaches in location theory, one which is purely deterministic (the median problem) and the 
other stochastic (the hypercube model). 


EXAMPLE 


The following example will illustrate some of our previous discussion. Suppose we want to 
locate three facilities on the simple network show.. in Fig. 1. The numbers next to the nodes are 
the fractions of demands from each node h,; j =1,...,5 and the numbers next to the links are 
the travel times. There are 6) possible distinct locations: 


{1,2,3}, {1,2,4},  {1,2,5}, {1,34}, (1,3, 5}, (1,4, 5} 
{2,3,4},  {2,3, 5}, {2,4, 5}, (3,4, 5}. 


The optimal location according to the standard 3-median problem is {1,2,5}, which can be 
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1 2 


Fig. 1. A simple 5-node network. 


obtained by hand. Suppose, however, that service requests occur in the network in a Poisson 
fashion with A = 4, and the service time for each one of the three units is exponential with 
identical means '=1. Let us assume that R=S units of time—the cost resulting when 
dispatching a reserve unit. We also assume that server preferences are determined solely by 
geographical proximity. 

The Jarvis algorithm with an initial location at the absolute 3-median, i.e. {1,2, 5}, converges 
after one iteration to the optimal solution at location {2,3,5}. The improvement achieved by 
moving from the location {1,2,5} to {2,3,5} is 3% in terms of the median problem with 
congestion. It is interesting to observe that the location {2,3, 5} is among the weakest possible 
locations in terms of the standard median problem (only location {1, 2, 3} is worse). This 
indicates that blind application of the absolute (deterministic) median problem to stochastic 
situations can lead to erroneous results, even for such simple networks. 
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